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1. Introduction 



Jet production at large transverse energies is a direct test of parton-parton scattering 
processes in hadron-hadron collisions. At large jet energy scales, the point-like nature 
of the partons can be probed down to distance scales of about 1CT 17 m by comparing 
data with QCD predictions. Within the experimental and theoretical uncertainties, 
data from the TEVATRON and CERN SppS generally show good agreement with the 
state-of-the-art theoretical next-to-leading order O (af) estimates based on massless 
parton-parton scattering over a wide range of jet energies [1, 2]. It is anticipated that 
the forthcoming Run II starting at the TEVATRON in 2001 will yield a dramatic 
improvement in the quality of the data with increased statistics and improved detec- 
tors, leading to a significant reduction in both the statistical and systematic errors. 
Subsequently, the start of data taking at the LHC will lead to a much enlarged range 
of jet energies being probed. 

It is a challenge to the physics community to improve the quality of the theo- 
retical predictions to a level that matches the improved experimental accuracy. This 
may be achieved by including the next-to-next-to-leading order QCD corrections 
which both reduces the renormalisation scale dependence and improves the match- 
ing of the parton-level theoretical jet algorithm with the hadron-level experimental 
jet algorithm. 

The full next-to-next-to-leading order prediction is a formidable task and re- 
quires a knowledge of the two-loop 2^2 matrix elements as well as the contribu- 
tions from the one-loop 2^3 and tree-level 2-^4 processes. In the interesting 
large-transverse-energy region, E? 3> ^ qU ark, the quark masses may be safely ne- 
glected and we therefore focus on the scattering of massless partons. For processes 
involving up, down and strange quarks, which together with processes involving glu- 
ons form the bulk of the cross section, this is certainly a reliable approximation. The 
contribution involving charm and bottom quarks is only a small part of the total 
since the parton densities for finding charm and bottom quarks inside the proton are 
relatively suppressed. We note that the existing next-to-leading order programs [1, 2] 
used to compare directly with the experimental jet data [3, 4] are based on mass- 
less parton-parton scattering. Helicity amplitudes for the one-loop 2^3 parton 
sub-processes have been computed in [5, 6, 7] while the amplitudes for the tree-level 
2^4 processes are available in [8, 9, 10, 11]. The parton-density functions are also 
needed to next-to-next-to-leading order accuracy. This requires knowledge of the 
three-loop splitting functions. At present, the even moments of the splitting func- 
tions are known for the flavour singlet and non-singlet structure functions F 2 and 
F L up to N = 12 while the odd moments up to N = 13 are known for F 3 [12, 13]. 
The numerically small Np non-singlet contribution is also known [14]. Van Neer- 
ven and Vogt have provided accurate parameterisations of the splitting functions in 
x-spa.ce [15, 16] which are now starting to be implemented in the global analyses [17]. 



1 



The calculation of the two-loop amplitudes for the 2 — > 2 scattering of light-like 
particles has proved more intractable due mainly to the difficulty of evaluating the 
planar and non-planar double box graphs. Recently, however, analytic expressions 
for these basic scalar integrals have been provided by Smirnov [18] and Tausk [19] as 
series in e = (4— D)/2. Algorithms for computing the associated tensor integrals have 
also been provided [20] and [21], so that generic two-loop massless 2^2 processes 
can in principle be expressed in terms of a basis set of known two-loop integrals. 
Bern, Dixon and Kosower [22] were the first to address such scattering processes and 
provided analytic expressions for the maximal-helicity-violating two-loop amplitude 
for gg — > gg. Subsequently, Bern, Dixon and Ghinculov [23] completed the two- 
loop calculation of physical 2 — > 2 scattering amplitudes for the QED processes 
e + e~ — > fi + fi~ and e + e" — > e"e + . 

In an earlier paper [24], we derived expressions for the two-loop contribution to 
unlike quark scattering, qq — > q'q', as well as the crossed and time reversed processes. 
The infrared pole structure agreed with that predicted by Catani [25] and we pro- 
vided explicit formulae for the finite parts in the s-, t- and w-channels in terms of 
logarithms and polylogarithms. Matrix elements for the other parton-parton scat- 
tering processes remain to be evaluated. In this paper we extend the work of [24] 
to describe the case of identical quark scattering. We use the MS renormalisation 
scheme and conventional dimensional regularisation where all external particles are 
treated in D dimensions to provide dimensionally regularised and renormalised ana- 
lytic expressions at the two-loop level for the scattering process 

qq -> qq, 

together with the time-reversed and crossed processes 

q + q -> q + q, 
q + q —> q + q. 

As in the unlike quark case, we present analytic expressions for the infrared pole 
structure, as well as explicit formulae for the finite remainder decomposed according 
to powers of the number of colours N and the number of light-quark flavours Np. 
For the contributions most subleading in N, there is an overlap with the two-loop 
contribution to Bhabha scattering described in [23] and the analytic expressions 
presented here provide a useful check of some of their results. 

Our paper is organised as follows. We first establish our notation in Section 2. 
The results are collected in Section 3 where we provide analytic expressions for the 
interference of the two-loop and tree-level amplitudes as series expansions in e. In 
Section 3.1 we adopt the notation of Catani [25] to isolate the infrared singularity 
structure of the two-loop amplitudes in the MS scheme. We give explicit formulae for 
the pole structure obtained by direct evaluation of the Feynman diagrams and show 
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that it agrees with the pole structure expected on general grounds. The finite O (e°) 
remainder of the two-loop graphs is the main result of our paper and expressions 
appropriate for the qq — > qq and qq — > qq (qq — > qq) scattering processes are are 
given in Section 3.2. Finally Section 4 contains a brief summary of our results. 



2. Notation 



For calculational convenience, we treat all particles as incoming so that 

q(Pi) + q(P2) + q(ps) + q(PA) -> 0, (2.1) 
where the light-like momentum assignments are in parentheses and satisfy 

pZ+pZ+pZ+p^O- 

As stated above, we work in conventional dimensional regularisation treating all 
external states in D dimensions. We renormalise in the MS scheme where the bare 
coupling a is related to the running coupling a s = a s (fi 2 ) at renormalisation scale 
(j, via 

po ^ + n% pi 



(2.2) 



«o S e = a s 

In this expression 

S e = (4vr) e e- e7 , 7 = 0.5772 . . . = Euler constant (2.3) 

is the typical phase-space volume factor in D = 4 — 2e dimensions, and (3 , Pi are the 
first two coefficients of the QCD beta function for N F (massless) quark flavours 

1 1C A - AT R N F n \1C\- 10C A T R N F - QC f TrN f 



Po 



Pi 



6 6 
For an SU (N) gauge theory (N is the number of colours) 



N 2 



C A = N, 



Tr = -t. 



2N J 1 ~~ 1 " 2 
The renormalised four point amplitude in the MS scheme is thus 

\M) = 47ra.[(|A<< >) - \M {0) }) + (g) - \M {1) ) 

\M^}-\M {2) })+0(al) 



(2.4) 



(2.5) 



2^ 



(2.6) 



where the \Ai^) represents the colour-space vector describing the i-loop amplitude 
for the s-channel graphs, and the t-channel contribution is obtained by ex- 

changing the roles of particles 2 and 4: 

\M (t) ) = \M (i) )(2^4). (2.7) 
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The dependence on both renormalisation scale /i and renormalisation scheme is im- 
plicit. 

We denote the squared amplitude summed over spins and colours by 

(M\M)=J2\M(q + q^q + q)\ 2 

= A(s, t, u) + A(t, s, u) + B(s, t, u), (2.8) 

where the Mandelstam variables are given by 

s = (pi+p 2 ) 2 , t = (p 2 + p 3 ) 2 , u =( Pl +p 3 ) 2 . (2.9) 

The squared matrix elements for the qq — ► qq process are obtained by exchanging s 
and u 

\M(q + q -> q + q) \ 2 = A(u, t, s) + A(t, u, s) + B(u, t, s). (2.10) 

The function A is related to the squared matrix elements for unlike quark scat- 
tering 

A(s,t,u) =J2\M(q + q^q' + q')\ 2 (2.11) 
A(t,s,u) =Y / \M(q + q'^q' + q)\ 2 (2.12) 

while B(s, t, u) represents the interference between s-channel and t-channel graphs 
that is only present for identical quark scattering. 

The function A can be expanded perturbatively to yield 



A(s,t,u) = 16n 2 a 2 s 



where 



A%s, t,u)+( ^) A%s, t, u) + (%L) 2 A 8 (s, t,u) + 



2tt7 v ' ' ' V2tt 



(2.13) 



't 2 + u 2 



A\s,t,u) = (M^\M (0) ) = 2(iV 2 - 1) (^— ^ - ej , (2.14) 

A e {s,t,u) = ((M {0) \M {1) ) + (M (1) |.M (0) >) , (2.15) 
A 8 (s,t,u) = ((M {1) \M {1) ) + (M (0) \M {2) ) + (M {2) \M {0) )) . (2.16) 

Expressions for A 6 are given in Ref. [26] using dimensional regularisation to iso- 
late the infrared and ultraviolet singularities. Analytical formulae for the two-loop 
contribution to .4 s , (M {0) \M {2) ) + (M {2) \M {0) ), are given in Ref. [24]. 
Similarly, the expansion of B can be written 



B(s,t,u) = 16ttV 



B%s, t, u) + (g) B\s, t, u) + (g) 2 B 8 (s, t,u) + 0( 



(2.17) 
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where, in terms of the amplitudes, we have 




(2.18) 



(2.19) 



B 8 (s,t,u) 




(2.20) 



As before, expressions for £> 6 which are valid in conventional dimensional regular- 
isation are given in Ref. [26]. Here, in order to complete the calculation of the 
two-loop contribution to quark-quark scattering, we concentrate on the next-to-next- 
to-leading order contribution B 8 and in particular the interference of the two-loop 
and tree graphs. 

As in Ref. [24], we use QGRAF [27] to produce the two-loop Feynman diagrams to 
construct either \M^) or | A^*' 2 '*). We then project by (M^\ or (M^ \ respectively 
and perform the summation over colours and spins. Finally, the trace over the Dirac 
matrices is carried out in D dimensions using conventional dimensional regularisa- 
tion. It is then straightforward to identify the scalar and tensor integrals present 
and replace them with combinations of the basis set of master integrals using the 
tensor reduction of two- loop integrals described in [20, 21, 29], based on integration- 
by-parts [30] and Lorentz invariance [31] identities. The final result is a combination 
of master integrals in D = 4 — 2e. The basis set we choose comprises 




(2.21) 




(2.22) 




(2.23) 




(2.24) 




(2.25) 



Pboxi(s, t) 



(s,t) 



(2.26) 




(2.27) 
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Xbox 2 (s,t) = { X ( s ^) ( 2 - 28 ) 

and 1 



Pbox 3 (s,t) = © (s,t), (2.29) 



where © represents the planar box integral with one irreducible numerator asso- 
ciated with the left loop. The expansion in e for all the non-trivial master integrals 
can be found in [18, 19, 20, 21, 28, 29, 33, 34, 35]. 



3. Results 

In this section, we give explicit formulae for the e-expansion of the two-loop contri- 
bution to the next-to-next-to-leading order term B 8 (s,t,u). To distinguish between 
the genuine two- loop contribution and the squared one- loop part, we decompose B 8 

as 

£8^8(2X0)^8(1X1)^ (31) 

The one-loop-square contribution B 8 ( - lxl ) is vital in determining B 8 but is relatively 
straightforward to obtain. For the remainder of this paper we concentrate on the 
technically more complicated two-loop contribution B 8 ( 2x0 ). 

As in Ref. [24], we divide the two-loop contributions into two classes: those that 
multiply poles in the dimensional regularisation parameter e and those that are finite 
as e — > 

B 8 {2x0) (s, t, u) = Voles + Finite. (3.2) 

Voles contains infrared singularities that will be analytically canceled by the in- 
frared singularities occurring in radiative processes of the same order (ultraviolet 
divergences are removed by renormalisation). 



1 Reference [20] describes the procedure for reducing the tensor integrals down to a basis involving 
the planar box integral 



Pbox 2 (s, t) = 



(s,t) 



where the blob on the middle propagator represents an additional power of that propagator, and 
provides a series expansion for Pbox 2 to O (e°). However, as was pointed out in [32], knowledge of 
Pboxi and Pbox2 to O (e°) is not sufficient to determine all tensor loop integrals to the same order. 
Series expansions for Pbox3 are relatively compact and straightforward to obtain and are detailed 
in [34, 35]. Pbox2 can therefore be eliminated in favor of Pbox3. We note that this choice is not 
unique. Bern ct al. [23] choose to use the Pboxi and Pbox 2 basis, but with the integrals evaluated 
in D = 6 — 2e dimensions where they are both infrared and ultraviolet finite. 
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3.1 Infrared Pole Structure 



Catani has made predictions for the singular infrared behaviour of two-loop ampli- 
tudes. Following the procedure advocated in [25], we find that the pole structure in 
the MS scheme can be written as 



Voles = -2 Re 



i(7W (0) |/ (1) (e)/ (1) (e)|-M (0) ) - - (M (0) \I {1 \e)\M^) 



+ (M {0) \I (1 \e)\M Wiin ) 

+e-^^l^ (7 + K ) <^ ( V (1) (2 e )|.M (0) > 
+ (M (0) \H^(e)\M {0) ) + {s <->*) 



where the constant K is 



K = (— - q C A - —TftNp. 
V 18 6 / 9 



(3.3) 



(3.4) 



In Eq. (3.3), the symmetrisation under s and t exchange represents the additional 
effect of the s-channel tree graph interfering with the t-channel two-loop graphs. 

The colour algebra is straightforward and we find that the s-t symmetric contri- 
butions proportional to 



N 



1 - e) - 



u 
~st 



(3.5) 



are given by 

(M m \I {1) (e)\M^} = (M {0) \M^} 

' 1 

TV 



on 



x 



1 3 
r(l-e) Ke 2 + Ye 



s 



2\< ! 

+ 



N 



2\ « 



N 2 + 1 
N 



.0. 

u 



2\ e 



(3.6) 



(M {0) \I {1 \e)I^(e)\M w ) = (M m \M (U) ) 



(0)\ 



r (0) i 



and 



2e 7 



r(i-e) 



1 3 \ 2 { N A — 3N 2 — 2 

^7 + 



+ 



3iV 2 + 1 
N 2 

1 

N 2 



2e 
u 

,2\ 2e 



2e 



AT2 

(iV 2 — 2)(iV 2 + 1 

2\ 2 M 



(7w (0, |ff( 2 )( e )|^W) = (7w (0) |^(°)) 

,2\ 2e 



=«7 



X 



2er(l-e) 



H (2) 



.E. 

s 



s 



t 



2c 



t 



u 



2c- 



(3.7) 



(3.f 
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where the constant is 



i 7 (i) + 3C F K + ^( 2 [3 C F - ^(3 C F 



/16 

IT 



7(^J C f Ca 



(3.9) 



Here ( n is the Riemann Zeta function, ( 2 — tt 2 /6, C3 = 1-202056 . . . and 

7(i) = (-3 + 24( 2 - 48C 3 ) C 2 F + - y C2 + 24( 3 ) C F C A + (J + ^2) CpTj^Np. 

(3.10) 

The square bracket in Eq. (3.8) is a guess simply motivated by summing over the 
antennae present in the quark-quark scattering process and on dimensional grounds. 
Different choices affect only the finite remainder. 

The bracket of 1^ between the t-channel tree graph and the finite part of the 
s-channel one-loop graphs is not symmetric under the exchange of s and t and is 
given by 



(A7 (0) |/ (1) (e)|^ (1)fin ) = 



3«7 



2\ e 



r(i-c) 

2\ e 



+ 



N 
'N 2 



+ 



N 



N 



.E 
s 



2\ e 



.EL 

' t 



N 



7 2 + 2~t 
N 2 + l 

2\* l 



t 



+ 



N 



.E. 

u 

A 

u 



2\ 6 



F\{s,t,u) 
(N 2 -l)F 2 (s,t,u) 



(3.11) 



The functions T\ and T 2 appearing in Eq. (3.11) are finite and are given by 
N 2 - 1 



Fi(s,t,u) 



2N 2 
1 

2e (3 - 2e) 

e e7 / 1 
r(l-e) le 



(N 2 -2)f 1 (s,t,u) + 2f 2 (s,t,u) 
N 2 — 1 

! (i - If - -h £ I - 

N 



+ 



N 
3 

2~e 



( 6 - 7e - 2e 2 ) - 1 (lOe 2 - 4e 3 ) 



Bub{s)(M i0) \M^) 



x 



N 



\ s J N 
1 3(1 -e) 
e ~ 3 - 2e 



u 



Bub(s) 



(7W (0) |^(°)) 



(0) l A/f(°)N 



(3.12) 



and 



F 2 (s,t,u) 



N 2 



1 r 



2iV 2 



fi(s,t,u) - f 2 (s,t,u) 



r(l -e) Ve 2 2e 



1 3 
^7 + ^ 



1 

iV 



1/, 



1 

N 



.EL 
t 



(M {0) \M^) 



(3.13) 
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with 



/i(s, t,u) = -£(l- 2e) [u 2 + t 2 - 2e (t 2 + s 2 ) + e 2 s 2 ] Box 6 (s, t) 
2u 2 - e (5s 2 + 6t 2 + 9s*) + (2s 2 + At 2 + st) e 5 



st 
2 



+ (s 2 + 3st) e 3 - ste 4 



Bub(s) - Bub(t) 



(3.14) 



2 

/ 2 (s, *,«) = -(!- 2e) [2m 2 - e (t 2 + s 2 + w 2 ) + 3e 2 s 2 + s 2 e 3 ] Box 6 (s, «) 



+ - 



st 



2m 2 - e (6s 2 + 6t 2 + lOst) + (3s 2 + 4t 2 + 3sf) e 2 
+ (s 2 + 2st) e 3 - ste 4 



Bub(s) - Bub(w) 



(3.15) 



These expressions are valid in all kinematic regions. However, to evaluate the pole 
structure in a particular region, they must be expanded as a series in e. We note 
that in Eq. (3.3) these functions are multiplied by poles in e and must therefore be 
expanded through to O (e 2 ). In the physical region u < 0, t < 0, Box 6 (w, t) has no 
imaginary part and is given by [23] 

\2 



Box 6 (w,t) = 



e e T (1 + e) T (1 - e) 
2sr (1 - 2ej (1 - 2ej 

+2e Li 3 (:r) - L x Li 2 (x) 



X L 3 — L 

3^ 2 ^ 

1 



1-4 l r 3 



-2e 2 I LU(x) + L,Li 3 (:r) - ^ 2 Li 2 (x) - ^L* x - - a LiL y + -L 2 X L 2 



7T 2 r2 7T 2 7T 4 \ 

"T x ~~ ~ 45 J + ^ ^ 



+ ( e 3 ) , 



(3.16) 



where x = —t/s, L x = log(x) and L y = log(l — x) and the polylogarithms Li n (z) are 
defined by 



^Li n _x(i) for n = 2, 3, 4 
o £ 



Li n (z) 

/•z //+ 

U 2 (z) = - -log(l-t) 

JO t 



(3.17) 
(3.18) 



Analytic continuation to other kinematic regions is obtained using the inversion 
formulae for the arguments of the polylogarithms (see for example [29]) when x > 1 

1 ^2 



Li 2 (a; + iO) = -Li 2 (-J - - log 2 (:r) + y + ivr log(x) 

(1 \ 1 7T 2 %1X 

-J - - log 3 (x) + y log(x) + y log 2 (a;) 

Li 4 (x + iO) = -LU (£) - ^ log 4 (x) + y log 2 ^) + ^ + y log 3 (x). (3.19) 
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Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 

The leading infrared singularity is O (1/e 4 ) and it is a very stringent check on 
the reliability of our calculation that the pole structure obtained by computing the 
Feynman diagrams agrees with that anticipated by Catani through to 0(l/e). We 
therefore construct the finite remainder by subtracting Eq. (3.3) from the full result. 

3.2 Finite contributions 

In this subsection, we give explicit expressions for the finite two-loop contribution to 
B 8 , Finite, which is given by 

Finite = -2Re ({M<® |AT (2)fin > + (M^^M^^ . (3.21) 

The identical-quark processes probed in high-energy hadron-hadron collisions are 
the mixed s- and t-channel process 

q + q -»• q + q, 

controlled by B(s, t, u) (as well as the distinct quark matrix elements A(s, t, u) and 
A(t, s,u) as indicated in Eq. (2.10)), and the mixed t- and -u-channel processes 

q + q ^ q + q, 
q + q -> q + q, 

which are determined by the B(t,s,u). We need to be able to evaluate the finite 
parts for each of these processes. Of course, the analytic expressions for different 
channels are related by crossing symmetry. However, the master crossed boxes have 
cuts in all three channels yielding complex parts in all physical regions. The analytic 
continuation is therefore rather involved and prone to error. We therefore choose 
to give expressions describing B 8 (s,t,u) and B 8 (t,s,u) which are directly valid in 
the physical region s > and u, t < 0, and are given in terms of logarithms and 
polylogarithms that have no imaginary parts. 

Using the standard polylogarithm identities [36] we retain the polylogarithms 
with arguments x, 1 — x and (x — l)/x, where 

t u x — 1 u , n 
x = —, y = — = l-x, = --. 3.22 

S S XI 

For convenience, we also introduce the following logarithms 

L x = log (—) , L y = log (—) , L s = log , (3.23) 



S J \ s 
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where \x is the renormalisation scale. The common choice fi 2 = s corresponds to 
setting L s = 0. 

For each channel, we choose to present our results by grouping terms according 
to the power of the number of colours iV and the number of light quarks N F , so that 
in channel c 

Timte c = 2 (^j^J (n 2 A c + B c + -^C c + NN F D C + ^E c + NpF^j . (3.24) 

Here c = st (ut) to denote the mixed s- and t-channel (u- and t-channel) processes 
respectively. 



3.2.1 The process qq — > qq 

We first give expressions for the mixed s-channel and t-channel annihilation process, 
qq — > qq. We find that 



A st = 



2LU(y) -2Lu(x) +2Li 4 



+ - 2 L x + 12 Li 3 (y) + 4 L y L x U 2 (y) 



( 23 \ /23 5 

+ - — -2L X + 4L„ Li 3 (a;)+ — + 12 L y + 2 L 2 + - vr 2 I Li 2 (a;) 
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1 



2 j ' ' ; . I 22 22 r ^ r | 11 ^ 2 | 22 r2 oo r : 592^ ^ 



x 3 
14 5 
T + 3 



27 



11 r2 31 13 a 



89 



695 \ 



+ -L? ; + 6L,j + -vr" + - L„ + — tt z -6C3 + ^ 



6 L ,+ 9 
347. 121 4 

C 3 7T 

18 S3 360 



9 / 

169 1 



216 / 



61 



18 6 



1673 \ 
108" y 



23213 
1296 



7T 



St 



+ 



4Li 4 (:r) + 24Li 3 (?/) + h L x + 12 \ U 3 (x) + ( - - vr 2 + 24 L y - 12 L x J Li 2 ( 



-? + ^)^+( -2(3-^^ + 12^ + 5^)^ 



12 



5 1 
2 + 3 



29 
If 



+-^-"-n 2 -12( 3 -4L y n 2 



u 

- + 

s 



3L 2 +3tt 2 + 3L 2 L y 



^7 + 



3L 



-32Li 4 (y) - 32Li 4 
-( -2Ll+( 



x — 1 



X 



+ 8L X Li 3 (y) + 2 - 28 L y + 18 L x Li 3 (x) 
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Some of these results overlap with the analytic expressions presented in Ref. [23] for 
the QED process e + e~ — > e + e~. To obtain the QED limit from a QCD calculation 
corresponds to setting Ca — 0, Cf — 1, T R — 1 as well as setting the cubic Casimir 
C 3 = (iV 2 — l)(iV 2 — 2)/N 2 to 0. This means that we can directly compare E st (cx 
CfTrNf) and F st (oc T^Np) but not C st which receives contributions from both C 3 
and C F . We see that (3.29) and (3.30) agree with Eqs. (2.50) and (2.51) of [23] 
respectively. 

The other coefficients, A st , B st , C st and D st are new results. 



3.2.2 The process q + q — > q + q 



The mixed t- and -u-channel process, q + q—>q + q is fixed by B 8 (t, s,u). We find 
that the finite two-loop contribution is given by Eq. (3.24) with 
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As in Section 3.2.1, we can compare some of these results with the analytic expres- 
sions presented in Ref. [23] for the QED process e + e~ — > e + e~, and we see that (3.35) 
and (3.36) agree with Eqs. (2.55) and (2.56) of [23] respectively. 

The other coefficients, A ut , B ut , C ut and D ut represent new results. 



4. Summary 

In this paper we discussed the two-loop QCD corrections to the scattering of two 
identical massless quarks. For the annihilation process, both s-channel and t-channel 
graphs are present. The interference of s-channel tree and two-loop graphs is deter- 
mined by v 4. 8( - 2x0 - ) (s, t, u) which is the same function that describes distinct quark scat- 
tering in the s-channel. Similarly, the interference of the t-channel tree and two-loop 
graphs is fixed by ^l 8< - 2x0 - ) (t, s, u). Explicit expressions for ^4 8 ( 2x0 ) are given in [24]. 
The modification to the matrix elements due to the interference of the s-channel 
tree graph with the t-channel two-loop graphs (and vice versa) is represented by 
# 8 ( 2x0 )(s, t, u) (see Eq. (2.10)). To obtain B s{ - 2x0 \ we have used conventional dimen- 
sional regularisation and the MS renormalisation scheme to compute the interference 
of the tree and two-loop graphs summed over spins and colours. 

The pole structure for £> 8 ( 2x0 ) (s, t, u) is given in Eq. (3.3) while expressions for the 
finite parts are given for the mixed s- and t-channels and mixed u- and t-channels in 
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Sees. 3.2.1 and 3.2.2 respectively. Together with the analogous expressions for unlike- 
quark scattering given in Ref. [24], they complete the analytic formulae required to 
describe the two-loop contribution to quark-quark scattering through to O (e°). 

These results form an important part of the next-to-next-to-leading order predic- 
tions for jet cross sections in hadron-hadron collisions. However, they are only a part 
of the whole and must be combined with the tree-level 2 — > 4, the one-loop 2 — > 3 as 
well as the square of the one-loop 2 — > 2 processes to yield physical cross sections. For 
the most part, the matrix elements themselves are available in the literature. Each 
of the contributions is divergent in the infrared limit and a systematic procedure 
for analytically canceling the infrared divergences needs to be established for semi- 
inclusive jet cross sections. Recent progresses in determining the singular limits of 
tree-level matrix elements when two particles are unresolved [37, 38] and the soft and 
collinear limits of one-loop amplitudes [39, 40], together with the analytic cancella- 
tion of the infrared singularities in the somewhat simpler case of e + e~ — > photon+jet 
at next-to-leading order [41], suggest that the technical problems will soon be solved 
for generic 2-^2 scattering processes. There are additional problems due to initial 
state radiation. However, the recent steps taken towards the determination of the 
three-loop splitting functions [12, 13, 14] are also promising. We therefore expect 
that the problem of the analytic cancellation of the infrared divergences will soon 
be addressed thereby enabling the construction of numerical programs to provide 
next-to-next-to-leading order QCD estimates of jet production in hadron collisions. 

Acknowledgements 

C.A. acknowledges the financial support of the Greek Government and M.E.T. ac- 
knowledges financial support from CONACyT and the CVCP. We gratefully ac- 
knowledge the support of the British Council and German Academic Exchange Ser- 
vice under ARC project 1050. This work was supported in part by the EU Fourth 
Framework Programme 'Training and Mobility of Researchers', Network 'Quantum 
Chromodynamics and the Deep Structure of Elementary Particles', contract FMRX- 
CT98-0194 (DG-12-MIHT). 

References 

[1] S.D. Ellis, Z. Kunzst and D.E. Soper, Phys. Rev. Lett. 69 (1992) 3615 [hep- 
ph/9208249]. 

[2] W.T. Giele, E.W.N. Glover and D.A. Kosower, Phys.Rev. D52 (1995) 1486 [hep- 
ph/9412338]. 

[3] F. Abe et al., CDF Collaboration, Phys. Rev. Lett. 77 (1996) 438 [hep-ex/9601008]. 
[4] B. Abbott et al., DO Collaboration, Phys. Rev. Lett. 82 (1999) 2451 [hep-ex/9807018]. 



18 



[5] Z. Bern, L. Dixon, D.A. Kosower, Phys. Rev. Lett 70 (1993) 2677 [hep-ph/9302280]. 

[6] Z. Bern, L. Dixon and D.A. Kosower, Nucl. Phys. B437 (1995) 259 [hep-ph/9409393]. 

[7] Z. Kunszt, A. Signer and Z. Trocsanyi, Phys. Lett B336 (1994) 529 [hep-ph/9405386]. 

[8] J.F. Gunion and J. Kalinowski, Phys. Rev. D34 (1986) 2119; 
S.J. Parke and T.R. Taylor, Nucl. Phys. B269 (1986) 410; 
F.A. Berends and W.T. Giele, Nucl. Phys. B294 (1987) 700; 
M. Mangano, S.J. Parke and Z. Xu, Nucl. Phys. B298 (1988) 653. 

[9] Z. Kunszt, Nucl. Phys. B271 (1986) 333; 

S.J. Parke and T.J. Taylor, Phys. Rev. D35 (1987) 313. 

[10] J.F. Gunion and Z. Kunszt, Phys. Lett. 159B (1985) 167. 

[11] J.F. Gunion and Z. Kunszt, Phys. Lett. 176B (1986) 163. 

[12] S.A. Larin, T. van Ritbergen, and J.A.M. Vermaseren, Nucl. Phys. B427 (1994) 41; 
S.A. Larin, P. Nogueira, T. van Ritbergen, and J.A.M. Vermaseren, Nucl. Phys. B492 
(1997) 338 [hep-ph/9605317]. 

[13] A. Retey and J.A.M. Vermaseren, hep-ph/0007294. 

[14] J. A. Gracey, Phys. Lett. B322 (1994) 141 [hep-ph/9401214]. 

[15] W.L. van Neerven and A. Vogt, Nucl. Phys. B568 (2000) 263 [hep-ph/9907472]; 
hep-ph/0006154. 

[16] W.L. van Neerven and A. Vogt, Phys. Lett. B490 (2000) 111 [hep-ph/0007362]. 

[17] A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne, hep-ph/0007099. 

[18] V.A. Smirnov, Phys. Lett. B460 (1999) 397 [hep-ph/9905323]. 

[19] J.B. Tausk, Phys. Lett. B469 (1999) 225 [hep-ph/9909506]. 

[20] V.A. Smirnov and O.L. Veretin, Nucl. Phys. B566 (2000) 469 [hep-ph/9907385]. 

[21] C. Anastasiou, T. Gehrmann, C. Oleari, E. Remiddi and J.B. Tausk, Nucl. Phys. 
B580 (2000) 577 [hep-ph/0003261]. 

[22] Z. Bern, L. Dixon and D.A. Kosower, JHEP 0001 (2000) 027 [hep-ph/0001001]. 

[23] Z. Bern, L. Dixon and A. Ghinculov, hep-ph/0010075. 

[24] C. Anastasiou, E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeomans, hep-ph/0010212. 

[25] S. Catani, Phys. Lett. B427 (1998) 161 [hep-ph/9802439]. 

[26] R.K. Ellis and J.C. Sexton, Nucl. Phys. B269 (1986) 445. 

[27] P. Nogueira, J.Comput.Phys. 105 (1993) 279. 



19 



[28] C. Anastasiou, E.W.N. Glover and C. Oleari, Nucl. Phys. B572 (2000) 307 [hep- 
ph/9907494]. 

[29] C. Anastasiou, E.W.N. Glover and C. Oleari, Nucl. Phys. B575 (2000) 416, Erratum- 
ibid. B585 (2000) 763 [hep-ph/9912251]. 

[30] K.G. Chetyrkin, A.L. Kataev and F.V. Tkachov, Nucl. Phys. B174 (1980) 345 ; 
K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159. 

[31] T. Gehrmann and E. Remiddi, Nucl. Phys. B580 (2000) 485 [hep-ph/9912329]. 

[32] E.W.N. Glover and M.E. Tejeda-Yeomans, Nucl. Phys. Proc. Suppl. 89 (2000) 196 
[hep-ph/0010031]. 

[33] R.J. Gonsalves, Phys. Rev. D28 (1983) 1542; 

G. Kramer and B. Lampe, J. Math. Phys. 28 (1987) 945. 

[34] T. Gehrmann and E. Remiddi, Nucl. Phys. Proc. Suppl. 89 (2000) 2512 [hep- 
ph/0005232]. 

[35] C. Anastasiou, J.B. Tausk and M.E. Tejeda-Yeomans, Nucl. Phys. Proc. Suppl. 89 
(2000) 262 [hep-ph/0005328]. 

[36] K.S. Kolbig, J.A. Mignaco and E. Remiddi, B.I.T. 10 (1970) 38. 

[37] JM. Campbell and E.W.N. Glover, Nucl. Phys. B527 (1998) 264 [hep-ph/9710255]; 
S. Catani and M. Grazzini, Phys. Lett. 446B (1999) 143 [hep-ph/9810389], Nucl. 
Phys. B570 (2000) 287 [hep-ph/9908523]. 

[38] F.A. Berends and W.T. Giele, Nucl. Phys. B313 (1989) 595; 

S. Catani, in Proceedings of the workshop on New Techniques for Calculating Higher 
Order QCD Corrections, report ETH-TH/93-01, Zurich (1992). 

[39] Z. Bern, V. Del Duca and C.R. Schmidt, Phys. Lett. 445B (1998) 168 [hep- 
ph/9810409]; 

Z. Bern, V. Del Duca, W.B. Kilgore and C.R. Schmidt, Phys. Rev. D60 (1999) 116001 
[hep-ph/9903516]; 

S. Catani and M. Grazzini, hep-ph/0007142. 

[40] Z. Bern, L. Dixon, D.C. Dunbar and D.A. Kosower, Nucl. Phys. B425 (1994) 217 
[hep-ph/9403226]; 

D.A. Kosower, Nucl. Phys. B552 (1999) 319; 

D.A. Kosower and P. Uwer, Nucl. Phys. B563 (1999) 477 [hep-ph/9903515]. 

[41] A. Gehrmann-De Ridder and E.W.N. Glover, Nucl. Phys. B517 (1998) 269 [hep- 
ph/9707224]. 



20 



